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Abstract 

The leading term for the energy of a bound state of fc-quarks and /c-antiquarks 
is proportional to its separation L. These fc-string configurations have a Liischer 
term associated with their quantum fluctuations which is typically a 1/L correction 
to the energy. We review the status of tensions and Liischer terms in the context 
of lattice gauge theory, Hamiltonian methods, and gauge/gravity correspondence. 
Furthermore we explore how different representations of the /c-string manifest them- 
selves in the gauge/gravity duality. We calculate the Liischer term for a strongly 
coupled SU(N) gauge theory in (2 + 1) dimensions using the gauge/gravity corre- 
spondence. Namely, we compute one-loop corrections to a probe D4-brane embed- 
ded in the Cvetic, Gibbons, Lii, and Pope supergravity background. We investigate 
quantum fluctuations of both the bosonic and the fermionic sectors. 



*christopher-doran@uiowa.edu 

tlpandoz@umich.edu 

* vincent-rodgers@uiowa.edu 

§kstiffle@gmail.com 



1 Introduction 



Since the emanation of Quantum Chromodynamics as the theory of the strong nuclear 
force, serious theoretical and predictive challenges appeared due to the strongly coupled 
nature of the theory at nuclear energies. What has emerged from these challenges is 
a technical and computational acumen that has advanced our understanding of other 
theories as well. Lattice gauge theory, advances in the Hamiltonian formulation of gauge 
theories, and string theories (including the gauge/gravity correspondence) are three ex- 
amples of these powerful theoretical tools that have emerged from this quest. At present 
the ability for theorists to match experimental data in the strongly coupled regime is 
very far from the success achieved in the perturbative regime. For this reason, theorist 
sometimes resort to comparing theories among themselves for states that are suitably 
accessible. For example the mass gap of the pure Yang Mills sector has been observed 
by the lattice community for some time and through string theory as well. For string 
theory, Polyakov gave a description of the mass gap by including the extrinsic curvature 
to Nambu-Goto action [I], where one can show that asymptotic freedom implies a mass 
gap. The question of the mass gap for 2 + 1 dimensions was also addressed using the 
Hamiltonian formulation in [2], [3]. Through this one sees that certain string theories and 
gauge theories belong to the same universality class. Suitable configurations that are cal- 
culationally accessible to all of these methods give theorists a standard for comparison. 
The fc-string epitomizes such configurations. In this note, we examine the gauge/gravity 
dual for a 2 + 1 dimensional fc-string configuration using the Cvetic, Gibbons, Lu, and 
Pope [I] (CGLP) supergravity background. Because this background is dual to a 2 + 1 di- 
mensional gauge theory, this configuration lends itself well to comparisons of both lattice 
gauge theories and Hamiltonian approaches to gauge theories. 

Here we give a brief overview of fc-strings and an outline of this note. For a more 
complete review of fc-strings, see j5]. fc-strings are a configuration of SU(N) color sources 
which result from k color sources in the fundamental representation stretched a large dis- 
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tance L from k anti-color sources. The gauge/gravity correspondence has been used to 
explore the relationship between certain configurations of low energy supergravity back- 
grounds and low energy ^-strings 0, [7j El El HQ] • Many of these supergravity calculations 
have been done for D3-branes, which are dual to four dimensional Yang-Mills theories. 
On the gauge theory and lattice gauge theory side of the correspondence, much of the 
focus has been in three dimensions [HJ [12l EU [HI [15], [T6|, [TTJ . Therefore the best test of 
the gauge/gravity correspondence through fc-string configurations is to find supergrav- 
ity backgrounds which are dual to 3d gauge theories. We will review in this paper the 
Cvetic, Gibbons, Lii, and Pope (CGLP) type IIA solution which is one such supergravity 
background [3]. 

From the gauge theory side of things, the classical energy of strongly coupled fc-strings 
is found to follow either a sine law or a Casimir law: 

Ek oc TfcL A;-string energy (1-1) 
T fc ociVsin— sine law (1.2) 

or 

Tfc oc k(N — k) Casimir law. (1.3) 
The precise form of the tension appears to be inconclusive [HI E2l [131 EU EES, ESI EH Ell 

Eg. 

In this paper, we will investigate specifically the gauge theory dual to iV stacked 
D2-branes, and stacked fractional D2-branes sourcing the CGLP type IIA supergravity 
background [1]. In this background, the low energy spectrum has been found to have a 
slightly lower tension than both the Casimir and sine laws [7]. 

We probe this calculation further, by calculating the one loop quantum correction to 
the low energy classical solution of the supergravity. Using vanishing boundary conditions 
on the probe D-brane, we find the energy correction is dual to an SU(N) Liischer term 
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with vanishing boundary conditions (i.e., the /c-string ends on immovable color sources), 

^Liischer ~~ ~~QL' (1-4) 

Our results constitute a prediction for the Liischer term for the fc-string. When comparing 
this result with the Liischer term coming from the fundamental string found in [201 12H 
[22], viz. — 7r/24L, one should be mindful to the fact that the /c-string configuration in 
the gauge/gravity correspondence has its origins in a large N setting where k/N is held 
fixed. In the same vein, naively setting k = 1 in the gauge/gravity correspondence may 
not recover the fundamental string. 

1.1 /c-String Ground State 

The fc-string tension provides the first level of comparison for various theories. In the 
context of M = 2 supersymmetric theories, Douglas and Shenker [23], examined the N- 
extended monopole condensation model of Seiberg and Witten for SU (2) [21] and found 
a spectrum of string tension that obeys a sine law, T k ~ iVsin(^). In a precursor to 
AdS/CFT, Hanany, Strassler and Zaffaroni were able to reproduce this spectrum of meson 
by using an M-theory fivebrane approach to QCD (MQCD) [25]. The lattice community 
has actively studied the fc-strings tension for some time [2S1 [221 12Z]- One of the issues 
that arises is whether these configurations exhibit a Casimir-like scaling 

T k « fc(l - k/N) 

for large N or a sine law where 

T k oc iVsin (^). 

Both of these behaviors respect the iV-ality but on the one hand, the 1/N expansion of 
QCD agrees with the sine law scaling [321 [IS], while on the other hand lattice calculations 
of [T7] favor Casimir scaling (1/N). 
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1.1.1 The Wealth of Information in 2 + 1 Dimensions 

Because the analysis is tractable, examining the fc-string configurations in 2+1 dimensions 
gives a wealth of information. In [2] Karabali, Kim and Nair predict that the string 
tensions in SU (N) gauge theories should be proportional to the quadratic Casimir of the 
representation of the flux. Based on the idea of effective dimensional reduction driven 
by a highly disordered vacuum, it was conjectured a long time ago [HI [12} [13] that this 
might also hold in D = 2 + 1 and D — 3 + 1. There is some additional evidence for this 
hypothesis from calculations of the potential between charges in various representations 
of SU{3). 

For a given k the smallest Casimir arises for the totally antisymmetric representation, 
and this should therefore provide the ground state k-string tension: 

o k _ k(N-k) 

V f ~ N-l ■ (i ' 5j 
This is the part of the Casimir Scaling hypothesis that we shall be mainly testing in 
this paper. For this purpose it is useful to have an alternative conjecture that possesses 
the correct general properties. A convenient and well-known example is provided by the 
trigonometric form 



(?k _ sin N 



/,'TT 

(1.6) 
Of sin F 

that was originally suggested on the basis of an M-theory approach to QCD [25J and has 
appeared in the context of the gauge/gravity correspondence. In fact the full prediction 
of Karabali et al. for o k is more specific than Eq.f ll.5l) . since it also predicts a value for 
Of in terms of g 2 , and including this gives: 

vr = e , (1.7) 

where Ca is the quadratic Casimir invariant for the adjoint representation defined as 
jacdjbcd = c A S ab . This may be written as: 

o k 1 k{N-k){N+l) 



(g 2 N) 2 8tt iV 2 
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Recent improvements to the Hamiltonian prediction of Karabali et al. [2B] have moved 
the value of oj to within —.3% to —2.8% of the lattice value. 

1.1.2 Representations for A;-strings 

In [29], [30], Gomis and Passerini studied various representations of Wilson loops in 
M = 4 supersymmetric Yang-Mills as related to particular representations of the gauge 
group. Inspired by this, we postulate that different supergravity backgrounds in the 
gauge/gravity correspondence can be identified with symmetric and antisymmetric rep- 
resentations of the /c-string. For the case of 2 + 1, both the Cvetic, Gibbons, Lii, and 
Pope (CGLP) [1] and the Maldacena-Nastase (M-Na) [31] supergravity backgrounds are 
relevant. We identify the following configurations with the group representations: 

• the probe D4 brane with world volume flux wrapping a three sphere in the CGLP 
background = antisymmetric representation; 

• the probe D3 brane with world volume flux wrapping a two sphere in the M-Na 
background = symmetric representation. 

We postulate that in the holographic context, at least in 2 + 1 dimensional theories, it 
seems that both types of formulas are possible: Casimir for CGLP and sine law for M-Na. 
Moreover, this suggests that the antisymmetric representations for the holographic case 
(CGLP) are "Casimir-law" while the symmetric representations (M-Na) are "sine-law". 
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Tk/Tf from Various Methods 
S=symmetric, A=antisymmetric, M=mixed, *=antisymmetric 



Group 


k 


CGLP 


MNa(Sine) 


Casimir 


lattice 


Karabali-Nair 


SU{4) 


2 


1.310 


1.414 


1.333 


1.353(A) 
2.139(S) 


1.332(A) 
2.400(S) 


SU(5) 


2 


1.466 


1.618 


1.5 


1.528* 


1.529* 




2 


1.562 


1.732 


1.6 


1.617(A) 


1.601(A) 




2.190(S) 


2.286(S) 


SU{6) 










1.808(A) 


1.800(A) 




3 


1.744 


2.0 


1.8 


3.721(S) 


3.859(S) 












2.710(M) 


2.830(M) 




2 


1.674 


1.848 


1.714 


1.752* 


1.741* 


SU{8) 


3 


2.060 


2.414 


2.143 


2.174* 


2.177* 




4 


2.194 


2.613 


2.286 


2.366* 


2.322* 



Table 1: Comparison of /c-string tensions from various methods. The values quoted 
are T^/Tf, where T& is the fc-string tension, and Tf is the fundamental string tension, 
i.e., k = 1. The CGLP tension is calculated from the transcendental Eqs. (12.54112.531) ; 
MNa(Sine) from Eq. (12. 551) ; Casimir from Eq. (11.31) . The A, S, and M data are calculated 
directly from [32], the * data is quoted directly from [17] . 

Table [T] shows this solution for supergravity backgrounds with lattice calculations by 
Bringholtz and Teper and Yang-Mills Hamiltonian calculations by Karabali et al. All cal- 
culations are in 2+1 dimensions where the quarks are in the anti-symmetric, symmetric, or 
mixed representations. Table [1] clearly shows that the CGLP tensions are lower than the 
Casimir law, lattice data and Karabali et al. calculations, and that the M-Na tensions are 
all higher. The full Table in 2 + 1, including all possible supergravity solutions with their 
corresponding brane embeddings, is not yet known. For example,other configurations to 



6 



be considered are the D6 branes in the CGLP background which would correspond to a 
symmetric state, D5 brane in the MNa background [31] (antisymmetric) . and the Witten 
2+1 QCD (nonextremal D3-brane compactified on a circle). As far as the explicit values 
of the tensions, our suggestion for which supergravity configuration/fc-string representa- 
tion correspondence is not well represented by the table. Both the CGLP and the M-Na 
tensions are more consistent with the antisymmetric representations. Further studies are 
needed to validate this issue and determine how these geometries manifest the /c-string 
tensor representations. The question still remains as to where the mixed representations 
would fit into the supergravity scenarios. 

1.2 Beyond the ground state: The Liischer term 

A detailed study of the flux tube between a quark and an anti-quark is an important 
window into the physics of confinement. A penetrating approach in this study is to 
consider the effective action for such string as an expansion in derivative terms. Liischer 
and Weisz considered all terms allowed by symmetries and built an effective theory for 
the excitations [33J and then went on to study the influence of the various terms on 
observables. Aharony and Katzburn were able to take this up to six derivative couplings 
and show that the full spectrum of the theory only depends on two free parameters 
[3l], the string tension a and a regularization dependent mass \i. The quark- ant iquark 
potential can be computed and shown to be 

V(L) = aL + f i + ^ + Oil/ L 2 ), T = -£(d - 2), (1.9) 

for large distances L, where d is the dimension of space-time. The 1/L correction in this 
formula is a quantum effect and can be used to determine a universality class for a large 
family of strings [21]. One remarkable property of the above expansion is that it sets in 
for relatively small values of L, for example it was argued in [33] that for values around 
0.5fm the above expansion is already valid. 
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Furthermore, there is an exponential reduction of the statistical errors in the lattice 
calculations of Liischer and Weisz [35] which allows them to compute the quark- ant iquark 
potential with high accuracy. Their results are completely consistent with theoretical 
model within a 0.04% and 0.12% depending on the quark-antiquark separation. 

We wish to compare the Lusher term in Eq. (jl.9p to that calculated from holographic 
models. The fundamental string is easy to find in most holographic models. It corre- 
sponds to a static fundamental string whose main contribution to the energy comes form 
the IR region of the dual supergravity background. Analyzing the massless modes of 
such fundamental strings in various confining backgrounds is rather simple since it is 
equivalent to counting symmetries. Since only those massless modes contribute to the 
Liischer term we conclude that: 



7 / d w = _|_ (CGLP H], M - Na [3T] backgrounds), 

fund ^ 



<ytj nd = - — ,(KS [36], M-Nunez [37] backgrounds). (1.10) 

As emphasized before, the above results come only from the massless modes. Namely 
in 3-dimensional models we have only one massless mode viz: X 1 , i = 2. Similarly, in 
4-dimensional models, there are two massless modes: X % with % — 2, 3, while all other 
modes are massive, including the fermions. 

Another very interesting configuration in confining theories is the one formed by k 
quarks separated a distance L by k anti-quarks. This configuration, known as the fc-string, 
also contains an intricate structure of excitations. We have studied such configurations 
at one-loop level in various supergravity backgrounds. We find that for both the KS and 
CGLP backgrounds, the four- dimensional*] and three-dimensional theories have the same 



Liischer term: 



k— string k— string ^ (1 11^ 



1 Here we correct a factor of a half missing in [ID] . 
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There are four massless modes for each. Let us explain: from the holographic point 
of view the difference comes from the fact that we have the standard massless Goldstone 
mode for the fluctuations in the transverse direction to the corresponding brane con- 
figuration, whereas the extra massless modes with respect to those of the fundamental 
string arise from the gauge field that is necessary to include holographically to get a 
representation characterized by k. In the case of the KS background our computation 
shows that there are two extra massless degrees of freedom coming from the fluctuations 
of the gauge field |10j . We verify later in this paper that in the three-dimensional case, 
the gauge field contributes three massless degrees of freedom to the Luscher term. 

2 /c-strings from type II Supergravity Duals 
2.1 The Basic Idea 

In this section, we give a brief sketch of how to calculate low energy /c-strings from a dual 
supergravity theory. The fc-string is thought of, from the supergravity dual, as a probe 
Dp-brane with electromagnetic charge in its world volume 

F = dA = F tx dtAdx + F H old A d<j>, (2.1) 

embedded in a supergravity background. 

This probe Dp-brane will wrap, or be tangent to, p + 1 out of the 10 bosonic super- 
gravity coordinates, X^. The remaining 9 — p supergravity coordinates will act as scalar 
fields, with dynamics on the probe Dp-brane. These scalar fields will enter into the action 
for the probe Dp-brane through the dilaton, $, and the pullbacks of the other bosonic 
supergravity sources 

Fn+i = dC n , H 3 = dB 2 . (2.2) 
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In addition, we will have fermionic fields with dynamics on our probe Dp-brane as 
well. The action we use for the probe Dp-brane, parametrized by (, is 

S Dp = -fi P I d p+1 ( e'V-det (g ab + F ab ) + ^ f e T aJ2 C i + Sf ( 2 -3) 

J J q 

F = B 2 + 2na'F, /i; 1 = (2tt)V (p+1)/2 . (2.4) 

where all of the fermionic fields are contained in Sf. The induced metric, g a b, is the 
pullback of the 10 dimensional supergravity metric 

_dX»dX» 

~~ d( a d( b 1 ' 

The generalization of this formula is how the n- forms listed in Eq. (12.21) are pulled back 
to the Dp-brane. 

For the low energy solution, we follow [7j and set the fermions all to zero, Sf = 0, 
consider bosonic fields constant on the Dp-brane coordinates, = Xq = constant, and 
integrate out the spatial Dp-brane coordinates from the action, Eq. (12.3[) . leaving us with 

S = J dt£(A a ,A a ,X£) (2.6) 

This particular solution leaves only the gauge field on the probe Dp-brane, A a with any 
dynamics with which to calculate the low energy Hamiltonian, which is proportional to 
the fc-string tension, T^: 

dC ■ 

H = -rj-A a — £ = LT k (2.7) 
2.2 Review of the CGLP Background 

Let us first review the CGLP type IIA supergravity solution. More details can be found 
in [H [7J [38]. The CGLP background is a solution with N coincident D2-branes 
In the string frame, the CGLP solution with these sources is found to be 

ds\ = G^dX^dX" = H- 1 ' 2 dx a dx p r) aP + H 1/2 ds 2 7 , (2.8) 
e" = gsH 1 ^ (2.9) 
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where rj a p is M 1,2 , and 



ds 2 = l 2 [h 2 dr 2 + a 2 (Drf) 2 + b 2 dtt 2 }, (2.10) 
X 2 = ^e ijk fi t Dfi i ADfi k , J 2 = ii\J\ X 3 = dX 2 = dJ 2 (2.11) 

In the above, I, m, and g s are constants, and a, b, h, Ui and H are functions of r. 

h 2 = (1-r- 4 )- 1 , G 2 = V(l-r- 4 ), 6 2 = V (2.12) 
Ul = r~ 4 + P(r)r~ 5 (r 4 - l)" 1 / 2 , w 2 = -^(r 4 - 1)^ + P(r)r^(r 4 - 1)~ 3 / 2 , 

-3 = rV-l)- 1 - 4r5(r4 _ 1)3/2 ^) (2-13) 

P M=f^= ( 2 - 14 ) 
m 2 Z" 00 

^ (r) = W J P( 2u *(P>^ ~ 3 Mp))dp. (2.15) 

The parameter I is similar to e in the deformed conifold [7J [39J, [36] and g s is the string 
coupling constant. 

The differential element Dfi 1 is 

Drf = dfj + e ijk A j fx k (2.16) 

where the fj, 1 are coordinates on a unitless M 3 constrained to a unit S" 2 surface, ///z 1 = 1. 
The fluxes associated with the solution are: 



H 3 = —a 2 mhdr A X 2 + —b 2 u 2 hdr A J 2 + —ab 2 u 3 X 3 , d = (2.17) 

F 4 = ^"Vx A c/PT 1 + m^7 1 G 4 , (2. 18) 

G 4 = a b 2 u 3 e ijk fj hdr A Drf A J k + a 2 b 2 u 2 X 2 A J 2 + \b^u x .] 2 A J 2 (2.19) 

where d 3 x = c?x° A c/x 1 A dx 2 . 

The A* are SU (2) Yang-Mills instanton one forms living on the S* 4 

A i = dQ% = (#, d X , dB, d(f>). (2.20) 
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and compose an anti-symmetric SU(2) Yang-Mills two form, J\ 

f = dA l + -e ijk A j A A k , (2.21) 
which satisfies the algebra of the unit quaternions, 

f P V P P = + ^kJip, (2-22) 

where g a p is the metric for an S" 4 : 

dn\ = gafsdn^d^ = dip + sin 2 ^{d X 2 + sin 2 X (d6 2 + sin 2 Odcj) 2 )). (2.23) 
As in [U [7], we select the J 1 to be 

J 1 = - sin ip dip A d% — sin 2 *0 sin 2 x sin 0d0 A d0 (2.24) 
J 2 = — sin?/; sin x A d6 — sin 2 ?/> sinx sin Odcp A g?x (2.25) 
J 3 = — sin?/; sinx sin #(£0 A d<p — sin 2 ?/> smxdx /\d9 (2.26) 

and we select the gauge where the solution for A % is 



A 1 = cos ipdx + cos i 
A 2 = cos ip sin xdO — cos x sin ( 
A 3 = cosip sinx sin 9d(p + cos x^ (2.27) 

We use the relations given in Eq. (12.21) to calculate B 2 and C3. In calculating B 2 , 
the identities in Eq. fl2.11l) are very helpful. Using these, we solve for B 2 , up to a total 
derivative, to be 



IB 2 = m 



£ hi^duj X 2 + m QT f 2 {r)duj J 2 



f^u) = a 2 (w)wi(tt)/i(w), f 2 (u) = b 2 (u)u 2 (r)h(u). (2.28) 

Notice that this vanishes when r = 1. This will be important for our ensuing calculations 
as it is where we will position our probe D4-brane. We choose the following solution for 

12 



C 3 : 

C 3 = _ x2drA df ^ dxl + ^eW^3 + ^u 2 (r)b(r) 2 a(r) 2 e ljk ^D^ A J k (2.29) 
g s tl[ry Hg s Zg s 

£(^>) = / sin 3 wdw, dfi 3 = sin 2 x sin 9 dx /\ d9 A d(p, (2.30) 

which is the same as that chosen in [7], up to a total derivative. 

The constant m is proportional to the number N of stacked fractional D2-branes 
that the background describes, which by the gauge/gravity correspondence, it is also the 
number N of colors for the dual supersymmetric SU(N) gauge theory. We can calculate 
the proportionality constant by using the Dirac quantization condition [7J 

F 4 = 87rV 3/2 iV, (2.31) 



's 4 

and the r — > 1 limiting behavior for G± 



G 4 ^-JdQ i . (2.32) 

Using these two conditions Eq.f l2.18l) gives m = 87ia' 3 ^ 2 g s N. 

Also, we must mention that there is another CGLP solution. To acquire the other 
solution, dVL A can be substituted for a metric over CP 2 . In this paper we will work only 
with the 4-sphere, dVL\. 

2.3 A Coordinate Transformation 

There is a singularity in the metric, Eq. (12.81) . at r = 1, because the function h(r), 
Eq.f 12.12]) . blows up here. This can be remedied by applying the coordinate transforma- 
tion 

T = yfi~T. (2.33) 
Now the point r = 1 is described by r = and the metric becomes 

ds 2 7 = I 2 (f(r) 2 dr 2 + a 2 (ZV) 2 + b 2 dn 2 ) 
f(r) = 2T (2.34) 

a/1 - (1 + T 2 )- 4 
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Notice that /(r) is finite as r — »• 0. Here are the expansions of all the aforementioned 
relevant functions in the r — > limit. 

a(r(r)) = r + 0(r 3 ), 6(r(r)) = -±=(1 + r 2 ), 

«i(r(r)) = | - 7r 2 + 0(r 4 ), u 2 (r(r)) = - \ + ^r 2 + 0(r 4 ), 
«3(r(r)) = ~ + ^r 2 + 0(r 4 ), /(r) = 1 + \r 2 + 0(r 4 ), 

777^ 777^ 7 

tf(r(r)) = tfo-# 2 r 2 + 0(r 4 )), tf = ^J , H 2 = —^ 

/oo 
p(2u 2 (p)u 3 (p) - 3u 3 (p))dp « 0.10693 . . . (2.35) 

One can use these expansions to show that _B 2 and C3 become, under the coordinate 
transformations Eq. fl2.33p . 

B 2 = --tJ 2 + 0(t 3 ) (2.36) 

of 

„ 2x 2 rdr A dx° A dx 1 3m ( ,. _ 1 9 _• , A ^. 

C3 = + 5- £ - ^T% k fi*D^ A J + O r 4 2.37 

g s H{r{T)Y 8g s \ 6 / 

2.4 Calculation of the /c-string tension 

We now outline the calculation of the k-string tension already presented by [7j. We label 
the world volume coordinates of the probe D4-brane as 

C = (t,x, X ,M), (2.38) 

where we are using the static gauge, and have fixed five of the bosonic supergravity 
coordinates to these D4-brane coordinates: 



X» = (t, x, x\ V, X, 0, 0, p 2 , /1 3 , r). (2.39) 

These 11 coordinates are really 10 independent bosonic coordinates, as the /i's are con- 
strained to (p*) 2 = 1. 
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We have then, the low energy action for the probe D4-brane, embedded in the CGLP 
background: 



-/14 J d 5 ( e'V-det (g ab + T ah ) + fi 4 JfAC 3 . (2.40) 



Here, B 2 and C3 have been pulled back to the D4-brane, the same way as the induced 
metric: 



dX» dX h 



9ab = QA „ a ^ G„ v . (2.41) 



We also turn on a U(l) gauge flux on the D4-brane 

27ia'Hl /2 F = E dtAdx (2.42) 

Since our D4-brane is placed at r = 0, where B 2 = 0, T becomes simply 

Hl /2 F = Hl /2 (2na'F) = E dt A dx (2.43) 

We examine the action Eq. (12.401) in static gauge, Eq. fl2.39[) . and at the classical solu- 
tion 

T = x 2 = /i 1 = /i 2 = /i 3 = 1 i) = xp . (2.44) 

Here we calculate the induced metric to be 

ds 2 Q = g a bd( a d( b = H~ 1/2 (-dt 2 + dx 2 ) + -dfi 2 (2.45) 

R 

whose scalar curvature is 

12 

At the classical solution, Eq. fl2.44p . Eqs. (12.36|2.37|) become 



R = esc 2 ^o- (2.46) 

Hi' I 2 



B 2 = 0, C 3 = cf = f^eWdOs- (2-47) 
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Plugging all of this into the action, Eq. fl2.4UI) . and integrating over x, x, 6, and 0, 
results in an effective action 

= [ dtC, (2.48) 



C = -aNL^l - E 2 sin 3 if) + qLNE^t/j) 

I 3 3 
a = — F , q = T-pra (2.49) 

2y/2irma' 2 3 / 2 / 1/2 
where L is the periodic length of the probe D4-branes x direction. Choosing a gauge where 
JF 01 = A x , leaves us with one conjugate variable with which to perform the Legendre 
transformation on the Lagrangian Eq. fl2.48p 

dC ■ 

H = -^A x - £ (2.50) 
8A X 

Because of the periodicity in the field E, the conjugate momentum to A x is quantized 
to an integer k [7j: 

2na'H 1 /2 ^ = ^- = kL. (2.51) 

BE 3A X V ' 

With this, we find the Hamiltonian, Eq. fl2.50p . to be 



H = aNL\l S m^ + ^(^-A , (2.52) 



Minimization with respect to ip results in the condition 

4fc w , , „af 



f faM + 3^ sin 2 cos (2.53) 
3N q z 

and the minimized Hamiltonian 

Hmin = TL = aNL sin 2 ^o^/sin 2 ip Q + (3a /q) 2 cos 2 ip 

a/q « 0.3083 (2.54) 

where T is the fc-string tension. Note that the parameter k is interpreted as the k 
quark- ant i- quark pairs in fc-strings. The tension, Eq. fl2.54p . and minimization condition, 
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Eg. (12. 531) . form a transcendental equation which can be solved numerically for given k 
and N. 

A similar, back of the envelopes calculation in the Maldacena-Nastase(M-Na) back- 
ground [31], leads one to a sine law: 

T~iVsin^. (2.55) 

Table [1] compares the tension calculated from the CGLP background Eq. fl2.54p to the 
sine-law, Eq. (12.551) . Casimir law, and various results from the lattice calculations and 
Hamiltonian formulation. 



3 Fluctuations of the Classical CGLP Solution 

Employing the same techniques as in [10], we will now fluctuate around the classical 
solution and calculate the one loop corrections to the classical energy, Eq. ( 12.541) . E±. 
Following [TUJ [3UJ |5T] , we calculate this correction to be given by the natural log of the 
functional determinate of the quadratic fluctuation of the partition function, Z 2 . 

e iElT = Z 2 = J DXDADQDQe iS2 (3.1) 

where S 2 is the part of the probe D-brane action quadratic in the fluctuations. 

To calculate this path integral, one would need to remove the gauge degrees of free- 
dom from symmetries such as U(l) gauge invariance of the gauge fields, diffeomorphism 
invariance of the probe D-brane, and re-symmetry of the supersymmetric D-brane action, 
via a Fadeev-Popov gauge fixing technique. We instead use the semi-classical techniques 
of [101 HH ED], and simply solve the gauge fixed equations of motion for the quadratic 
fluctuations, and sum over the resulting eigenvalues: 

Ex = EP + E[ f) (3.2) 
2 We thank A. Armoni for a discussion of this point and its overall relevance. 
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where E\ is the bosonic energy correction and E[ is the fermionic energy correction. 

We now seek to the find the equations of motion for the quadratic fluctuations. The 
general fluctuations will be given by 

X» = + 5X» bosons (3.4) 
A a = + 5A a U(l) gauge fields (3.5) 
6 = + 50 fermions (3.6) 

We will first concentrate on the bosonic fluctuations. For fermionic fluctuations, we 
will later look to Martucci and collaborators [4*2l l4*3l HH H5] to find the proper type IIA 
supergravity fermionic action. 

3.1 Bosonic Fluctuations 

Of the 10 independent bosonic coordinates, we fluctuate only five of them: the five not 
statically set to the D4-branes world volume: 

x 2 (() = + 6x 2 ((), ^(C) = ^o + #(C), r(C)=r + 5r(C), A0=ti + ¥(0 

(3.7) 

where the subscript zero, as in /Xq, refers to the classical value for the field, specified in 
Eq. (12.441) . Under this fluctuation, the induced metric becomes, to quadratic order in the 
fluctuations 

ds 2 = dsl + ds l + ds l (3-8) 
I 2 

dsl = H l' 2 2 sin ( 2 ^o)^3<¥ ( 3 - 9 ) 
.2 _ f rr-i/idSx 2 dSx 2 i /2 2 / d5r 85t IdStpdStpW b 



+ 



dS2 = ( ^ + H « l [WW + 2WW 1 1 dCdC + 



+ 5r 2 



J^^dt 2 + dx 2 ) + Hl ,2 l 2 {A^dn^ + sin 2 ^ (1 - ^")^3 



I 2 

+ (^ 2 # 1/2 - cos(2^ )^3 i = 1,2, (3.10) 
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B 2 becomes 

B 2 = B { 2 1] + B { 2 2) (3.11) 
B { 2 1] = — sin 2 sin X St d\ A d9 (3.12) 

B^ = —St ( 5ji l sin 2 ^ sin 2 X sin 9 d9 A d<ft + 5fi 2 sin 2 ipo sin x sin 9 d<p A d\+ 



dSip \ 

+ smip sinx sin 9 — — d( a A d<p + Sip sin(2'0o) sin% d\ A dd 1 , (3.13) 

and C3 becomes 

C 3 = Cf + cf+Cf , (3.14) 

Cf = -^ sin 3 Vo^ 3 , (3.15) 
°9s 

Cf > = f 9 sin 2 cos ^o^ 2 - ^ + ^ 5tA dQ 3 (3.16) 
16^ \ sin x sino 1 / 

The fluctuation, Eg. (13.41) . leads to a simple expansion for the U(l) gauge field: 

dt A dx + SF, 



2na'Hl /2 



SF = (d a 5A b )d( a A d( b . (3.17) 

As the dilaton depends on r through Eq. fl2.9p . the dilaton value, to second order in 
the tau fluctuations, becomes 

e -* = e -*o( 1 + ^ (Jr 2) ( 3 . 18 ) 
Ho 

e^=g s H l Q /A (3.19) 
Using Eq. (13.81) through Eq. fl3.18l) . we calculate the bosonic action, Eq. fl2.40p . to second 



19 



order in the fluctuations: 



S{ b) = J ^-det(g(eff))d 5 ( k d 
4 6 ) = _ y yC det(( 7 ( e //))c/ 5 c|c :E V a 5x 2 V a (5x 2 + 

+ c T [V a <5rV a 5r + m 2 ( X , #)5r 2 ] + c A 

+ tot a lder 1 vat 1 ves■ 



(3.20) 
(3.21) 



V a 6if>V a 6i(> - 



+ 



167T 



+ 



(3.22) 



where the covariant derivatives are with respect to g^ e ^\ an effective metric on the 
D4-brane 



1 

9xx 



ds 2 = 9 { ; b ff) dcdc b = —{-de + dx 2 ) + ^dni 



12 3 A 2 / e 



9a 



H 2 T 2 R 3 m 4 ' 



(3.23) 



where i? is the same scalar curvature as in Eq J2.46l The U(l) gauge current j a , and 
m T (x, 0) are 

V^sin 98t) 



J 



sinx sin 6 



sin 2 x 



ml {x, 9) = m ro + "? csc2 X esc 2 9, 
6 

and the various constants are 

fx^R 3 ^ 2 l 3 H min 2/ 4 



(3.24) 
(3.25) 



I 4 



3R „ 9HL jTlW 1 ' 2 



1" 



A8y/3g 8 A h™ 2 
Qi> = 



I m 2 T 32/ 1 /27r 3 ma /2 
I 9 

mm 



4/ 1/2 7r 2 12 2 a 



8v / 6A 3 7T 2 Z 15 /2l2 



5 a> 



I 



m 



t0 



16m/, 



3/2 



:i + 8/ 



7 
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i?, ci 



.4 



iA 9/2 / 2 m 4 



12 7 2tt 2 Z 3 



(3.26) 



Clearly, } is an integral over total derivatives. This confirms that we are truly 
fluctuating around a classical solution. The second order action, S% , as it stands, is 
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difficult to acquire eigenvalues from, as can be seen in the equations of motion from this 
action which are listed in Appendix [K\ There we solve one of the equations of motion, 
and discuss solutions of the others. 

To calculate the Luscher term, we integrate out the spherical degrees of freedom, x, 
9, and 0, which will leave us with the same number of massless modes as before. Since 
the Luscher term only depends on massless modes this process should lead us to the 
same Luscher term as would be calculated from the full five dimensional theory. This 
was explicitly found to be the case in [10] where the massless modes were independent of 
the angular degrees of freedom. 

To proceed with this integration, we consider the fluctuations to be independent of 
the S 3 variables, 



5X» = SX"{t, x), SA a = 5A a (t, x) 



(3.27) 



and we integrate out the S* 3 from the action Eq. fl3.22p . This results in an effective action 



R 



S$> ff = -V 3 J dt dx^c x V m 5x 2 V m 5x 2 + 
+ c T [V m 5TV m 5r + m 2 Te 5T 2 ] + c A 



V m c^V m c^ - ^g x M 2 



5F mn 5F mn + g xxJ m 8A r 



ca 



+ 77r(V m 6A x V m 6A x + 2W m 5A e W m 5A e + / 1 V m <5^V m 5^)+ 

107T 



+ total derivatives > , 



(3.28) 



where the indices m and n now sum only over the coordinates t and x, and are raised 
and lowered by the two dimensional Minkowski metric 



Vmnd( m d( n = -dt 2 + dx' 



(3.29) 



and the effective 5t mass, m Te , the constant V3, and the integral, I\ are 

3/2 



h 



esc 8 dO. 



(3.30) 
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The equations of motion of the action, Eg. (13.281) . are 



V z 5x z = 

V 2 ^ + - 9*M + ^9±SF tx = 



'ip 

V 2 5t - m 2 Te 5r = 
V 2 5A = 0, i = X ,0,<t> 
V m 5F mn = Airg xx3 n 

To solve these equations, we move to Fourier space 

(oo 2 - p 2 )5x 2 = 

(u; 2 - P 2 + - z^f^M, + P A t ) = 

(lu 2 -p 2 - m 2 Te )5r = 

(^ 2 - P 2 )^ = o, % = x^A 

p 2 SA t + puA x = -iAng^QippSi/j 
to 2 5A x + puA t = -i4ng xx Q,puj5ip, 

and work in temporal gauge, 5A t = 0. This leaves us with six eigenvalues 

p 2 4 fold degenerate 
p 2 + m 2 e , 
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(3.31) 
(3.32) 

(3.33) 
(3.34) 
(3.35) 



(3.36) 
(3.37) 

(3.38) 
(3.39) 
(3.40) 
(3.41) 



p 2 + m 2 



where 



R 
~2 



ca 
2c^ 



4TiQ 2 ,g : 



2 „2 



I 



45.9992— < ml < oo. 

m v 



(3.42) 
(3.43) 
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S (f) = 0± J d 5 Ce"*V-det MQT' D4 



6, (3.45) 



The calculation for the one loop correction to the bosonic fc-string energy, is given 
in Appendix [Bj the result for large quark separation L is 

Ef ] = -^ L - l -{m Te + m i ). (3.44) 

This is in contrast to the Liischer term found in [20], EH E2], which is — 7r/24L. 



3.2 Fermionic Fluctuations 

Following the series of papers [42j [43j HH H5], we use the K-symmetry fixed fermionic 
action for the probe D4-brane found in [45] : 

(M- l ) ab Y a Df ] + (M~ l ) ab Y h W a - A 

where the definitions are found in Appendix ICl 

Applying the fluctuation Eq. fl3.4p to the fermions leads us to an action second order 
in the fermionic fluctuations, 

(M~ l ) ab T a Df ] + (M- l ) ab T b W a - A 5Q. (3.46) 

As we did in the bosonic case, we integrate out the S 3 , resulting in the effective quadratic 
fermionic action 

S 2eff « J dtdx5QT' Di {{M- 1 ) mn T m d n + M f )6Q, m,n = t,x. (3.47) 

The matrix Mf is given in Appendix [C] We solve the Euler equation from this Lagrangian 
by Fourier transform 

r^ 4 (z(M- 1 r n r m p n + M / )5e, m,n = t,x pt = -u, Px = p. (3.48) 

The eigenvalue solutions to this equation, outlined in Appendix O, are 



S!P = yJ d 5 (e-*V-detM8er' D4 



LU = ±yp 2 + Oil ± «2 

±\/p 2 + a 3 ± a 4 



to 



±V a 7(P) + «S(P) ± <*6 (P) /0 , n x 

w = s > ( 3 - 49 ) 

±y/a 7 (p) -a 5 (p) ±a 6 (p) 
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These eigenvalues prove difficult to regulate when one calculates the one loop fermionic 
energy Eq. fl3.3p . This calculation has been postponed until a later time. 



4 Conclusion 

We analyzed the ground state and one loop quantum corrections for SU (N) fc-strings. For 
the ground state, we compared string theory results to lattice gauge theory and Yang-Mills 
theory results. We investigated the hypothesis that D4-branes in the CGLP background 
would describe quarks in the anti-symmetric representation and that D3-branes in the M- 
Na background would describe quarks in the symmetric representation. For the tensions, 
we found, in fact, that both of these were more closely related to anti-symmetric quarks. 
We concluded that more research must be done with other probes like D2-branes and 
D6 in the CGLP background and D5-branes in the M-Na background to find a clearer 
correlation with symmetric, anti-symmetric, or mixed quark representations. 

Furthermore, we analyzed the fluctuations from the ground state for SU(N) /c-strings 
in 2 + 1 dimensions from a duality relation with a D4-brane probing the type IIA CGLP 
supergravity background. We found equations of motion for bosonic and fermionic fluc- 
tuations that are relevant for the Luscher term. Interestingly, we found the Luscher term 
was — 7r/6L, the same as in our previous investigation for 3+1 fc-stringg^. 

We predict the Luscher term for the fc-string through the gauge/gravity correspon- 
dence. Since this result is based on large N with k/N fixed one can only make a qualitative 
comparison to the work of [201 [2U [22] • In fact all that we might be showing in this case 
is that the quantum excitations of the /c-strings are different from the lattice result for 
the fundamental string. 



! Hcre we correct a factor of a half missing in [10] 
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A Five Dimensional Bosonic Equations of Motion for 
the Fluctuations 

Applying the variational principle to the action, Eq. fl3.22p . results in the field equations: 

V 2 5x 2 = (A.l) 
V 2 # + f «ty + °-^SF tx = (A.2) 

V V - m 2 T ( X , 9)5t + S^L csc x SFg x = (A.3) 

V a 5F ab - 4vr/ = 0, (A.4) 

where V a is the covariant derivative compatible with Eq l3.23l 
The solution to Eq. flA.lj) is 

Sx 2 = J duodp ^ n ' l ' m) {p,u)e^ x -^Y nlm { X ,e^), (A.5) 

n>l>\m\ 

where the Y nlm (x, 0, <j)) are the spherical harmonics on an S 3 j45| 

Y nlm ( X ,6,<f>) = c nl ^^P^\ / / 2 2 (cos X )Y^(e, ( j ) ), 

/ (n + l)(n + / + lj! 
Cnl = \] {^1)1 ' (A ' 6) 

and P l n {x) are the associated Legendre polynomials. The S 3 spherical harmonics Y m (x, 

satisfy the eigenvalue problem 

V 2 Y nlm ( X , 6, <P) = -n(n + 2)Y nlm ( X , 6, 0), (A.7) 
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where V 2 is the Laplacian for an S 3 whose action on scalar functions such as Y nlm (x, 9, 4>) 
is explicitly given by 

^ 2 = *k xdx) + ^ e(8in m + a??"*) ■ (A - 8) 



With that said, solving Eq.l lA.ll) with the solution Eq.( 1A.5l) results in the eigenvalue 
problem 



— (to 2 -p 2 ) - ^n{n + 2) 
9 xx 6 



x = (A.9) 



The rest of the equations prove quite difficult and require perturbation theory to solve; 
their solution is not given here. 

B Calculation of the One Loop Bosonic Energy 

Following [TQl HQl Ej , we define the one loop bosonic energy as a sum over the bosonic 
frequencies, Eq. (13.421) . 

E P = \ E ( 4 P + vV + <, + y/p 2 + m 2 ) . (B.l) 

p 

To describe quark- ant i- quark sources affixed to the end of the string, we use vanishing 
boundary conditions 

p = rni/L (B.2) 
As in [TU], we use (- function regularization to regularize the massless modes 

oo 
n=l 

For the massive modes, instead of using the techniques of [17], we regulate the sum 
by using a cutoff 

oo [M] 

v 7 ™ 2 + m 2 — > J2 ^ n2 + M2 ( B - 4 ) 

n=l n=l 
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where [M] is the largest integer less than M. Now we are able to use the binomial 
expansion for the square root 

n \ 2 i 

m) 



E^^= M EEu ( 

n=l n=l g=0 

oo /ix [M] 

CO / 1 \ 00 

(7=0 n=l 



/I 



2q 



n 



2<J 



m E(^) m " 29 ^(°) 



M 

T 



(B.5) 



where we have again used ^-function regularization, and the approximation [M] — > 00 is 
valid for large M oc L, large L corresponding to large quark separation. 

Using these results, we find for large quark separation L, the one loop bosonic energy 

is 



dm 



C Fermionic Dp-brane Action 

Throughout this paper, and unless otherwise noted, Latin indices a,b,c, . . . are D-brane 
indices, with the exception of i, j, and k, which will always refer to the unit S 2 coordinates, 
/x\ Greek indices, a, f3,/i, u, . . . are 10 dimensional curved indices, and underlined Greek 
indices, /i, z/, a, . . . are 10 dimensional flat indices. 

The fermionic action for the quadratic fluctuations was found in [J2j H3], 0U 05] . It is 
given by 



S (f) = t± j d 5 (e-*V-detM5er' D4 



(Af 



-i)«*r o D< 0) + (M~ 1 ) ab r b w a 



A 



5e. asm 
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where 50 is a 32 component spinor, constrained by 

T U 5Q = 5Q. (C.l) 

where T 11 = r 0123456789 , and P^ 2 "^" is the totally antisymmetric product of gamma 
matrices. The flat gamma matrices satisfy a Clifford algebra 

{ra, r^} = 2^ (c.2) 

where i]^- is the 10 dimensional Minkowski metric. A useful consequence of this anti- 
commutation relation is the following identity 

rVa/3 = fj, ^ v ± a ^ /3, (C.3) 

which can be generalized to any number of gamma matrices. 

We define A = A« + A (2) , M = g + T, F 4 = F 4 + H 3 A C u F 2 = dC u where 

DT = d a + j^riV + j^H a ^ (C.4) 



(C.5) 



A« = i(r^* + ^ir^F"*) (C.6) 



a (2) = ^ - i^r^) 



(C.7) 



= _ V detff g) (rll)P /2+l V- ( ggg raxa,...^^ ^ (ag) 
e a 1 a 2 ...a p+1 

= ( P +l)\^d^ Taia2 - ap+1 (C,9) 

and g a i a 2--- a P+i i s a density, i.e., takes values of ±1, or 0. 

The 10 dimensional curved T^'s are related to the 10 dimensional flat T^'s by the frame 
fields, T M = e-^r^, and the D-brane T a 's are pulled back from the curved 10 dimensional 

rvr = 
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The frame-fields, 



G^ = eVi%i (CIO) 
for the CGLP background can be written in a 10 dimensional representation as 

e% = e\ = e\ = H^'\ e% = IfH^, 

e~ 3 = esc 9e~ 4 = esc x esc 9e~ 5 = esc ij) esc x csc 9e% — IbH 1 ^ 4 , 

e\ = lH^a%^ k ^\ e% = lH^acsc9^Aie^^, a = 4,5,6 (Gil) 
09 d(p 

with parametrization of the unit S 2 , (/x 1 ) 2 = 1, given by 

/i 1 = sin6 l cos0, /i 2 = sin6 l sin0, /i 3 = cos6 l (C.12) 
In the above, the the 10 independent bosonic coordinates are numbered ... 9 as 

X» = (t,x,x 2 ,ij, X ,9, (f>,9,4>,r) (C.13) 

We use the frame fields to calculate the spin connection for the 10- dimensional bosonic 
space 

%~ = %~ = rP&a (W + e £ T V) (C. 14) 

where the Christoffels, T a , and the inverse frame-fields, e M , are given by 

e£ = VE± G^e\ (C.15) 

With this, we calculate the spin connection for the CGLP background. Pulling its lowered 
index back to the D4-brane and evaluating it at the classical solution, we find the only 
non-vanishing components to be 

f2 x — = f2 x — = cos^o, — — costpo sinx, 

54 n 87 63 87 , • • Q 

\l e — = \l d — = cosx, = = cos^o sinx smu, 

fir- = fL— = cosx sin#, fir- = fL— = cos 9, (C.16) 
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With this we calculate the term in the action which contains the spin connection to be 

- (M" 1 )' 1 '' T a n b ^T^ = M C + cot X Mi + esc X cot 9 M 2 

1 [~R 

m c = -\/- cos^ (3r 3 + r 498 + r 579 + r 687 ) 

2 V o ~ 

M 1 = lw^(2r 4 + r 587 + r 697 ) 

2 V o ~ 

M 2 = ~W^(r 5 + r 698 ) (C.i7) 

2 V o ~ 

This is the only term in the action that has 9, x dependence, modulo the measure. Many 
of the formulas in Eq. (IC.4l) -Eq. (IC.9l) simplify 

W a = -le^^F^r^Ta (C.18) 
A (1) = j { H a ^\ A< 2 > = __i_e*°F Q ^r Q ^ (C.19) 



r' =i- ^_!^ r n (i - -r n r a ^ afc ), (c.20) 



all of which are, again, x, 9 independent. 

As in the bosonic case, we investigate S 3 independent solutions for the quadratic 
fluctuations 

8e = se(t,x) (c.21) 

leaving us with an action of the form 

S£ ff <x J dtdx J dxd9d(f)sm 2 x sin 95QT' D4 ((M~ 1 ) mn T m d n + Mf+ 

+ cot x Mi + esc x cot 9 M 2 )SQ, m,n = t,x, (C.22) 

where 

M f = M c + {M- l ) ab Qr a ^r^ + Y h WA - A (C.23) 

Integrating out the S 3 as before, it is easy to see that the terms proportional to Mi 
and M2 integrate to zero, leaving us with 

S 2eff oc J dtdx5QT' Di ((M- 1 ) mn T m d n + M f )6Q, m,n = t,x. (I3T47D 
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We solve the Euler equation from this action by Fourier transform 



Y'{i{M- l ) mn T mPn + M f )5Q, m,n = t,x Pt = -u, p x = p. QSM) 



We now pick a representation for the 32 x 32 gamma matrices 



T a = ia 1 <g> cr 1 <g> cr 1 ® a 3 ® a , r 1 = a 1 ® a 1 <g> a 2 ® a ® a 
r 2 - = a 1 ® a 1 g> a 1 ® a 1 <g> a 3 , r a = a 1 ® a 1 g> a 1 ® a 1 <g> a 2 



r 4 - = a 1 <g> cr 1 ® a 3 <g> cr° <g> a . 



= -a 1 <g> a 1 ® a 1 ® cr 1 <X) 



cr 2 ® ct^ ® ct^ ® ct^ (x) a . 



o- 1 <g> a 2 <g) a ® a 



o 



a 



T 1 = cr 1 <g> cr 3 ® a ® cr Q$ a 
r a = a 1 g> cr 1 <g> a 1 ® a 2 <g> a 4 



o 



(C.24) 



where ® means tensor product, and the <t m are the Pauli spin matrices, augmented with 
the identity: 



a 



1 
1 



a 



1 

1 



cr 



-i 

1 



a 1 



1 
-1 



(C.25) 



This representation leaves us with a diagonal T 1 !, and so through the constraint, 
Eq.f lC.lj) . we are able to set the lower 16 components of 5Q to zero. At the same time, 
this reduces the 32 equations in Eq. ( 13.481) to 16 independent equations. These equations 
can be reorganized into the following form 



u5Q = H f 5Q 



(C.26) 



where the Hamiltonian has the block diagonal form 

( H ± \ 
H f = H 2 



(C.27) 



V H 3 J 
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and Hi and H 2 are 4x4 matrices, and H 3 is an 8 x 8 matrix 



/ 



V 

where the c's are constants. 
The eigenvalues of Hf are 



p 







Ca 







p 


c 6 




,H 2 = 




c c 


-p 







Cd 







-P j 




~P 







-C c 








-p 


-c rf 


C 3 





Ci 




p 





c n 


-c b 


Q 





p 











-Cfc 





-P 











-Cfc 





Cn 








— Q — c e 





c n 








-c.f 



/ 



w = ±\/p 2 + ai ±«2 



V 



P Ci C e 
p Cf Ci 

Cj Cg ~P 

c h Cj -p 

\ 



-Cfc 

-Cfc 



-c n 



"P -C h ~Cj 





P 



±Vp 2 + « 3 ± a 4 

±\/«7(p) + a 5 (p) ± a^"(p) 
±a/ot(p) -« 5 (p) ±«e (p) 



(C.2* 



(C.29) 



± 



where ai, a 2 , a^, and are constants combinations of the c's in IC.28t and a^,a 6 , 
and a 7 are functions of p: 
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a 7 (p) 


= 4 0) + 4V, 




= E a?y + 4A(p), 

n=0,2,4 




= 2a§(p)-3# L (p)±A;(p), 


A(p) 


= 12 UW +A(P) > 


2& 3 (p) 


= P*(p) + y/ 0l(p) - 40l(p) 


ft(p) 


= E ^ V, 

n=0,2„8 




= E /^V, 

n=0,2„12 




n=0,2„6 



Here, the a- and /3\ n> are constant combinations of the c's from IC.28I 



?(«) 



(C.30) 
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